Non-saturating magnetoresistance of inhomogeneous conductors: comparison of 

experiment and simulation 
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The silver chalcogenides provide a striking example of the benefits of imperfection. Nanothreads 
of excess silver cause distortions in the current flow that yield a linear and non-saturating transverse 
magnetoresistance (MR). Associated with the large and positive MR is a negative longitudinal MR. 
The longitudinal MR only occurs in the three-dimensional limit and thereby permits the determina- 
tion of a characteristic length scale set by the spatial inhomogeneity. We find that this fundamental 
inhomogeneity length can be as large as ten microns. Systematic measurements of the diagonal 
and off-diagonal components of the resistivity tensor in various sample geometries show clear evi- 
dence of the distorted current paths posited in theoretical simulations. We use a random resistor 
network model to fit the linear MR, and expand it from two to three dimensions to depict current 
distortions in the third (thickness) dimension. When compared directly to experiments on Agjj-^Se 
and Ag2_|_^Te, in magnetic fields up to 55 T, the model identifies conductivity fluctuations due to 
macroscopic inhomogeneities as the underlying physical mechanism. It also accounts reasonably 
quantitatively for the various components of the resistivity tensor observed in the experiments. 

PACS numbers: 72.20.My,72.15.Gd,72.80. Jc 
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I. INTRODUCTION 

Macroscopic inhomogeneities were generally associated 
with a diminution of galvanomagnetic effects, due to the 
reduction of carrier mobility caused by impurity scatter- 
ing. Classical effective media theoriesi'^'^i^i^iii^ have 
shown, on the other hand, that conductors normally 
exhibiting no change in resistance R when a magnetic 
field H is applied, can possess a positive magnetore- 
sistance, AR/R EE (RiH) - i?(0))/i?(0), in the pres- 
ence of spatial fluctuations of the resistivity, provided 
that the characteristic length scale of the inhomogene- 
ity is much larger than the mean free path. This in- 
duced magnetoresistance depends generically on whether 
the magnetic field is oriented perpendicular (transverse 
MR) or parallel (longitudinal MR) to the current direc- 
tion. In recent experiments, Solin and co-workers2iiS 
have harnessed this effect by crafting composite materi- 
als in which current deflections around highly conducting 
inhomogeneities significantly enhance the intrinsic trans- 
verse MR of Hgj_^Cda;Te. This designer approach de- 
pends on the mismatch in the physical MR of the two 
components and leads to a significant, quadratic magne- 
toresistance that saturates hy H ^ 0.5 T. 

In addition to geometric enhancements to the galvo- 
magnetic response of materials with intrinsic physical 
MR, there is the possibility of designing materials and 
devices where the MR arises solely from the macroscopic 
inhomogeneities. In such a case, where the disorder- 
induced component of the magnetoresistance can be sep- 
arated from the physical magnetoresistance, the func- 
tional form of the MR not only can be tuned to be linear 
with magnetic field, but it is possible to craft a response 



that continues to at least a MegaGauss. Essential de- 
sign considerations involve the pertinent inhomogeneity 
length scale, which can involve both the dopant distri- 
bution and polycrystallinc structure, adjacency to a zero 
band gap condition to amplify the effects of fluctuations 
in the carrier mobilities, and the ease of growth and fab- 
rication. 

The silver chalcogenides stand out as favorable ma- 
terials for the investigation of large, linear and non- 
saturating MR because of their simplicity and tun- 
ability. Perfectly stoichiometric AgjSe and Ag2Te are 
non-magnetic, narrow-gap semiconductors whose elec- 
tron and hole bands cross at liquid nitrogen temper- 
atures. They exhibit negligible physical magnetoresis- 
tance>ii, as predicted from conventional theories. By con- 
trast, minute amounts of excess Ag or Te/Sc — at levels 
as small as I part in 10,000 — lead to a huge and lin- 
ear magnetoresistance over a broad temperature range, 
with no sign of saturation up to 60 T^^d^d'^d^i-'^^d'^ . In 
particular, the unusual linearity extends deep into the 
low field regime with H <^ I/Mi where ^ is the typical 
mobility of the material. Such behavior shows no resem- 
blance to conventional semiconductors, where the mag- 
netoresistance grows quadratically with field and reaches 
saturation at fields typically of order 1 T. Therefore, it 
has been argued that the observed magnetoresistance 
must be caused by the inhomogeneous distribution of 
excess/deficient silver ion o^^'^^ . The spatial inhomo- 
geneities in the silver chalcogenides result from the clus- 
tering of Ag ions at lattice defects or grain boundaries 
when quenched from the high temperature a-phase to the 
semiconducting /3-phase with lower Ag solubility. They 
may take the form of highly conducting nanothreads or 
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lamellae along the grain boundaries of the polycrystalline 
material's. 

Abrikosov was the first to point out the essential role 
played by inhomogeneities in the magnetoresistance ob- 
served in the silver chalcogenides, emphasizing the im- 
portance of the zero band gap condition. The phe- 
nomenon is quantum in origin and relies on a combi- 
nation of a linear dispersion and a single, partially-filled 
Landau bandi^i^i. This picture also has been applied 
to the large, positive MR in epitaxial Bi filmsS^ and in 
suitably doped InSb^. By contrast. Parish and Little- 
wood (PL) ^^i^^ demonstrated that purely classical, geo- 
metric effects of the inhomogeneities could be responsible 
for the non-saturating magnetoresistance in Agj^^-Se and 
Ag2-|-5Te. They introduced a two-dimensional random re- 
sistor network model capable of simulating macroscopic 
media with complicated boundaries. Their simulations 
of strongly inhomogeneous semiconductors derive a lin- 
ear MR from the Hall resistance picked up from macro- 
scopically distorted current paths caused by variations in 
the stoichiometry. Moreover, when both positive (hole- 
like) and negative (electron-like) values of the mobility 
are sampled, as will occur at bank crossing, the func- 
tional form of the MR is most linear (consistent with 
experimentsii) . The magnetic field scale for the onset 
of the linear response is then set by the width of the 
mobility distribution. These results also agree qualita- 
tively with subsequent studies of two-component media 
involving exact results and effective medium approxima- 
tions2^i2Sii^, which show that the non-saturating trans- 
verse magnetoresistance will be maximal when the high 
field Hall resistance crosses zero. 

The exact distribution of inhomogeneities is difficult to 
resolve via either scattering or imaging techniques for the 
relevant doping levels in the silver chalcogenides, but the 
longitudinal MR has been found to serve as an effective 
probe of the spatial extent of the excess silver nanoth- 
reads^^. The negative longitudinal MR and, in partic- 
ular, the emergence of an unexpected cross-term in the 
resistivity tensor coupling the longitudinal and transverse 
response, appears as a physical manifestation of distorted 
current paths, with an attendant inhomogeneity length 
scale. Missing at present, however, is the theoretical de- 
piction of current distortions in the thickness dimension, 
and the abnormal longitudinal MR associated with it. 
The majority of theoretical investigations on classical in- 
homogeneous conductors focus solely on the transverse 
MR, the major exception being the study of isotropic 
media with insulating inclusions, where the longitudinal 
MR is substantially larger than the transverse MRiiS. 
Moreover, the modeling has been restricted to transport 
properties in the bulk, whereas the presence of bound- 
aries is clearly important for the unusual longitudinal 
MR of the silver chalcogenides. 

In this article, we extend the PL random resistor model 
to three dimensions in order to permit an explicit com- 
parison to the longitudinal MR results. We find that 
simple realizations of macroscopic inhomogeneity repre- 



sent a crucial ingredient for understanding the full mag- 
netoresistive response observed in the experiments. The 
PL model allows one to simulate the various components 
of the resistivity tensor, and generate visualizations of 
the current distortion in the thickness dimension, results 
only accessible through the computer simulations. More- 
over, we test the theory and fix its parameters via direct 
comparison to the temperature and magnetic field de- 
pendence of the transverse MR. 

The paper is organized as follows. Section H de- 
scribes the sample preparation and measurement tech- 
niques, while Sec. HI reviews both the 2D network model 
and the 3D expansion we use to visualize the distorted 
current paths in the presence of a longitudinal magnetic 
field. In Sec. IV, we compare the simulation results 
and experimental data regarding the linear transverse 
magnetoresistance, the transverse-longitudinal coupling 
Pxz{H) and the negative longitudinal magnetoresistance. 
We present our conclusions in Sec. V. 



II. EXPERIMENTAL METHODS 

Appropriately weighted amounts of high purity Ag and 
Te (99.999%, Alfa Aesar) were melted to created poly- 
crystalline Ag2-j_5Te at desired stoichiometries, both sil- 
ver rich (n-type, 5 > 0) and silver deficient (p-type, 
(5 < 0). The compound was ground and loaded into out- 
gassed, fused silica ampoules inside a helium glove box, 
then baked at 50 K above the melting point (1170 K) to 
ensure complete mixing. Slowly cooled specimens with 
typical dimensions (4.0 x 1.0 x 0.4) mm^ were cut per- 
pendicular to the cylindrical axis to avoid possible dopant 
variations due to temperature gradients inside the fur- 
nace. Ag2-|-5Se polycrystals were prepared similarly from 
stoichiometric AgjSe, and appropriate weights of Ag or 
Se were added to reach the desired compositions. Sample 
thicknesses from 400 /im to 15 /xm were obtained by me- 
chanically shaving the same specimen with precise depth 
control under the microscope. 

We performed four-probe magnetotransport measure- 
ments using a conventional ac bridge technique in the 
ohmic and frequency-independent limits. Ultrasonically 
soldered InBi contacts were placed on the top, bottom 
and sides of the specimens to determine the pxx , Pzz and 
Pxz components of the resistivity tensor. The transverse 
and longitudinal MR was determined by averaging values 
at positive and negative field directions. The longitudi- 
nal MR is much smaller than the transverse MR and thus 
requires attention to alignment of the magnetic field and 
current directions (within ±5°). We used a 14/16 T su- 
perconducting magnet in a helium-3 cryostat for low field 
measurements and the 60 T short-pulsed magnet at the 
National High Magnetic Field Laboratory at Los Alamos 
for the high field response. The latter provides 16 ms 
long field pulses from the discharge of a 1 MJ capaci- 
tor bank and, when combined with the NHMFL-designed 
synchronous digital lock-in amplifier with response times 
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much shorter than commercial lock-ins. can provide low 
noise data during each 16 ms shot. 



III. THEORETICAL MODELS 

For classical magnetotransport, where the mean free 
path of the charge carriers is much less than the typical 
length scale of the disorder, we can take Ohm's law to be 
valid locally in space: 

E(r) = p(r) j(r). (1) 

Here, j is the current density, E is the electric field and 
p is the resistivity tensor. For simplicity, we will as- 
sume that the inhomogeneities are isotropic and that 
each point in space is characterized by a single charge 
carrier, so that the resistivity tensor acquires the follow- 
ing form in a magnetic field H = Hz: 

/ 1 (3 0\ 
p = pi -(3 1 ] (2) 
V 1/ 

Clearly, this tensor is dependent on just two parameters: 
the carrier mobility p (since the dimensionless variable 
/3 = pH) and the scalar resistivity p. A realistic inho- 
mogeneous semiconductor will generally possess multiple 
charge carriers, but their effect on p is typically small 
with Ap/p < l^a. 

In order to address the problem of strong inhomo- 
geneities, PL previously introduced a two-dimensional 
square network of four-terminal resistor unit a^^'^'^ , which 
we will briefly review here. As shown in Fig. 1(a), the re- 
sistor unit is taken to be a homogeneous disk whose trans- 
port properties can be calculated easily using Eq. (1). 
Four terminals are required to take account of the Hall 
voltage that develops perpendicular to the current di- 
rection in a transverse magnetic field. The voltage dif- 
ferences Vi between terminals are linearly related to the 
currents at each terminal via an impedance matrix , 
i.e. Vi = Zij Lj . Like Eq. (2), the impedance matrix is 
completely determined by the mobility p and an effec- 
tive scalar resistivity p. We can thus generate a square 
N X N random network (Fig. 1(b)) by connecting up the 
disks using Kirchoff's laws and then randomly varying p 
and p for each disk. Specifically, we consider the distri- 
bution of p to be Gaussian with average (p) and width 
A/z, while the distribution for p must obviously satisfy 
p > (e.g. we can take p = rf' ^ where 77 has a Gaussian 
distribution). Note that negative and positive values of 
p correspond to electron-like and hole-like carriers, re- 
spectively, and that the case where {p) = corresponds 
to an equal number of electrons and holes since our cho- 
sen distribution is symmetric about {p) =0. In order 
to determine the network resistance R, we apply a po- 
tential difference V across the network in Fig. 1(b) and 
then sum up the currents along the left (or right) edge, 
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FIG. 1: Two-dimensional resistor network model and its ex- 
tension to three dimensions. The original 2D resistor unit (a) 
is assumed to be a homogeneous disk with four equally spaced 
terminals, each of which is associated with a current ii and 
a voltage difference Vi . The resultant N x N resistor net- 
work (b) has a constant potential difference V applied from 
left to right, and no current can enter or exit the network at 
the top and bottom edges. Unlike the previous resistor net- 
work modeU^, we assume periodic boundary conditions for 
the currents on the left and right edges l/' = In the 

three-dimensional version, the basic unit is naturally a six- 
terminal element (c) of unit volume, but it now represents 
a direct discretization of Eq. (1) rather than any particular 
geometry such as a sphere. 



/ = SO that R = V/I. We avoid boundary ef- 

fects associated with perfectly conducting electrodes on 
the left and right edges of the network^! by taking peri- 
odic boundary conditions for the currents, = I^. 

In the limit iV — > 00, this model successfully simulates 
the transverse MR of an inhomogeneous semiconductor. 
However, in order to model the longitudinal MR, we must 
extend this network to three dimensions, even in the case 
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FIG. 2: Simple realizations of two-dimensional disorder consisting of two different resistivities pi and p2, where the magnetic 
field H is parallel to the applied voltage V . Patterns (a) and (b) have no longitudinal magnetoresistance, while (c) has a small, 
saturating one. The magnetoresistance of (d) can be tuned to be large by setting pi/ p2 2> 1. 



where the inhomogeneity is purely two-dimensional. We 
should emphasize that there are numerous subtleties as- 
sociated with adding a third dimension. First, the lo- 
cal resistivity tensor (2) is obviously not isotropic like 
it is in 2D owing to the fact that there is now a pre- 
ferred direction along the magnetic field. Indeed, this 
property is responsible for the standard phenomenon of 
current jetting--. In addition, the resistivity in 3D con- 
tains the dimensions of length (i.e. R = PxT") ^'^ 
is not scale invariant as it is in 2D. From a theoretical 
standpoint, the extra dimension also makes it impossible 
to use duality relation o^'^'^^i'^° or conformal mappings 
to obtain exact results in 3D, although much progress 
has been made with effective medium approximations for 
two-component systems^i^i^i^^, and by exploiting an anal- 
ogy with advection-diffusion problems in the continuous, 
weakly disordered cas o^i^i^^ . 

The above considerations mean that there is no 
straightforward generalization of the circular resistor unit 
to 3D. In fact, deceptively trivial geometries like a sphere 
or a cylinder with its axis perpendicular to the magnetic 
field have a non-saturating magnetoresistance^^. Instead, 
we must employ a "brute-force" approach by directly dis- 
cretizing Eq. (1) into six-terminal elements as in Ref . [ssl 
and then constructing 3D networks that simulate basic 
macroscopic inhomogeneities. Referring to Fig. 1(c), we 
set a voltage and current for each terminal such that each 
voltage is defined with respect to the element center and 
incoming currents are defined as positive. The impedance 
matrix relating voltages and currents is then given by: 

Z^, = ^{S,, + (3) 



where 

/ -1 1 0\ 
1 0-1000 




Vo 000/ 

A major advantage of this model over effective medium 
approximations is that we can include sample boundaries, 
a crucial ingredient for understanding the observed neg- 
ative longitudinal MR in the silver chalcogenides^^. In 
particular, we can simulate a four-probe measurement, 
where the voltage probes placed on top of the sample are 
separate from the electrodes at the ends of the sample 
that supply the current (by contrast to the two-probe 
measurements illustrated in Fig. 2). 

The spirit of our investigation will be to construct 
simple realizations of disorder that capture the salient 
features of experiment, thereby gaining insight into an 
otherwise difficult theoretical problem. We depict in 
Fig. 2 simple examples of possible two-dimensional in- 
homogeneities with two different resistivities pi and p2- 
Note that there is no problem with using perfectly con- 
ducting electrodes on the ends of the sample, since un- 
like in a transverse field, there is no contact MR in a 
longitudinal field. We see straightaway that the classi- 
cal longitudinal MR of Fig. 2(a) and (b) is zero since 
all the current is parallel to the potential difference V 
(or, equivalently, the magnetic field H) and is thus in- 
sensitive to the Lorentz force. However, once there is 
a non-zero current perpendicular to the magnetic field. 
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Hall voltages will build up in the z-direction, which, in 
turn, can affect the current flow. The chessboard pat- 
tern in Fig. 2(c) exhibits a small magnetoresistance for 
Pi 7^ P2 that eventually saturates once the current has 
found a path across the sample that is parallel to H . The 
ideal pattern is one where the current is forced to flow 
perpendicular to H and strongly experiences the Lorentz 
force as in Fig. 2(d) for pi/ P2 ^ 1- One may even re- 
gard this "percolation" pattern as representing a generic 
filamentary current path analogous to those seen in the 
2D network simulation a^^i^^ , although in this case the 
current is filamentary even at zero magnetic field. Sur- 
prisingly, differences in mobility /i between the two com- 
ponents have hardly any effect on these simple patterns, 
so we take /i to be uniform. 

For the remainder of this article, we take Fig. 2(d) 
to be our basic toy model for longitudinal magnetore- 
sistance. We have computed the magnetoresistance of 
different x Ly x network sizes and we find that 
it has converged for > 12, Ly > 16 and > 24 
networks considered throughout. 



IV. RESULTS 

A. Linear Magnetoresistance in a Transverse 
Magnetic Field 

In our initial studies of the silver chalcogenides we con- 
centrated on the galvanomagnetic behavior in a trans- 
verse magnetic field. Semiclassically, the dependence of 
resistance on transverse field is contained in the prod- 
uct UJT, where uj is the frequency at which the mag- 
netic field causes the charge carriers to sweep across 
the Fermi surface. Given a universal mean free time 
r oc 1/np, where n is the carrier density, the field de- 
pendence of the transverse MR can be collapsed onto a 
universal curve by expressing AR/R solely as a function 
of the magnetic field, carrier density and the resistivity: 
AR/R = fiiJT) = f{H/np). 

This functional dependence is known as Kohler's rule, 
obeyed by a great number of metals and semiconduc- 
tors, with a positive, quadratic MR expected to saturate 
for CUT « 1. By contrast, in doped silver chalcogenides, 
the anomalous transport is characterized by a large, lin- 
ear magnetoresistance that shows no sign of saturation 
at 60 T and extends down into the low field regime 
{H <C l//i) where one generally expects AR/R oc H^. 
Moreover, the curves at different temperatures do not 
coincide with the scaling form AR/R = /{ujt), but are 
shifted by a multiplicative factor that depends on tem- 
perature^. 

The PL model is able to account for the linearity at 
low fields by showing that the linear, transverse MR of 
a heavily inhomogeneous semiconductor crosses over to 
quadratic behavior at a field set by the inverse of the 
mobility distribution width Afi, rather than the average 
mobility (fi), provided that A/i/(/i} > 1 . In addition. 
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FIG. 3: Normalized transverse magnetoresistance, AR/R, as 
a function of magnetic field measured at T = 110 K. The solid 
line is a theoretical fit from 20 x 20 random resistor networks 
averaged over 10 realizations with (p) = and A^~^ = 0.58 
T, consistent with the PL hypothesis Afj,/{fi) > 1. Inset; As 
seen in a blow up of the low field region, the inverse of the 
mobility distribution width, Afj,~^ , sets the approximate scale 
for the crossover to conventional quadratic behavior. 



the character of the mobility distribution determines the 
magnitude of the linear response: at sufficiently large 
fields, AR/R oc ApH for strongly disordered semicon- 
ductors with A/i/ (fi) > 1 . Therefore, the galvanomag- 
netic behavior of the silver chalcogenides naturally devi- 
ates from a Kohler's plot, since its temperature depen- 
dence is no longer simply tuned by the change in carrier 
mobilities, but also by the disposition of the sample in- 
homogeneity. 

For a quantitative understanding of the temperature 
and field dependence of the transverse MR in silver 
chalcogenides, it is necessary to compare directly the the- 
oretical predictions of the PL model with experiment. 
Fits to the data will yield an estimate of the mobility dis- 
tribution in the silver chalcogenides and will test whether 
the crossover from quadratic to linear MR is set by the 
absolute value of the mobility fi or by the width of the 
mobility distribution A/i. We plot in Fig. 3 the trans- 
verse MR in n-type Ag2+5Se {S = 0.0001) at both low 
and high fields, fit to the random resistor network model. 
We find that (/i) — yields the best fit, with A/i~^ ~ 
0.58±0.05 T, confirming the expectation that the silver 
chalcogenides are in the regime of A/x/(/x) :s> 1, with a 
crossover field set approximately by Ap,~^. While the 
value for A/i~^ is still slightly higher than the observed 
crossover field, it is remarkably good given that we are 
considering a simple random resistor network model with 
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FIG. 4: Visualizations of currents and voltages for the two-dimensional "percolation" pattern depicted in (a), taking pij pi = 
10,000 and using slightly smaller networks of 12 x 16 x 24 elements for better visibility. The currents correspond to black 
arrows where the magnitude is given by the arrow size, while the voltages are represented by different colors, assuming the 
voltage drop across the sample is IV. Panel (b) shows the currents and voltages as viewed from the top in the absence of a 
magnetic field, where the transport is independent of z. In a large magnetic field /3 = 10, the current develops a component 
in the z-direction, resulting in a loss of current at the top surface as depicted in (c). Panel (d) shows a slice sX x = Lx/2 that 
reveals the behavior of the current in the z-direction. 



only two adjustable parameters. Our toy model is over- 
simplified in many ways: it is two-dimensional instead of 
three, it is restricted to a narrow range of disorder types, 
and it neglects the contact effect between network disks. 
Of particular interest is the fact that the electrons and 
holes are present in equal proportions, as indicated by 
(/i) = 0, consistent with the claim that the magnetore- 
sistance is most linear when both positive and negative 
values of the mobility can be sampled. For an inhomoge- 
neous system, however, the averaged mobility cannot be 
simply obtained from a measurement of the Hall resis- 
tance. The experimental determination of the mobility 
distribution requires a more precise characterization of 
the internal structure, as a real material generally con- 
tains a small fraction of heavily doped regions imbedded 
into a stoichiometric background with intrinsic carrier 
concentration. Indeed, an unambiguous theoretical pre- 
diction for the Hall resistance has proved elusive for the 
random resistor network in this regime^!. 



B. Current Distortions in a Longitudinal Field 

The MR in conventional semiconductors is caused by 
the curvature of the carrier trajectories in a magnetic 
field. In a longitudinal geometry, where the magnetic 
field lies in the same direction (the x-direction) as the ap- 
plied voltage, the carriers' momentum in the x-direction 
is conserved for a uniform current flow, yielding a nearly 
negligible longitudinal MR. The diagonal tensor compo- 
nents Rxz{H), representing transverse- longitudinal cou- 
plings, should be zero due to the vector nature of the 
magnetic field. When strong inhomogeneity exists, how- 
ever, current distortions in the y and z directions can 
significantly affect the galvanomagnetic measurements. 

While experiments cannot access local current pat- 
terns directly, fortunately one can numerically simulate 
them. Our simple two-dimensional percolation pattern 
in Fig. 4(a) provides an ideal starting point for demon- 
strating the existence of field-induced current distortions 
along the sample thickness. At zero field the transport 
is independent of z, due to symmetry, and there is no 
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FIG. 5: (a) The voltage drop measured across the thickness 
of n-type Ag2Se, reveahng distorted current flows as a direct 
manifestation of sample inhomogeneities (depicted in the in- 
set). A conventional six-probe technique was used to compare 
signals from the edge (E) and middle (M) of the specimen, (b) 
Numerically simulated Rxz{H) components of a 18 x 24 x 
30 three-dimensional network (illustrated in the inset). The 
potential on each side of the sample was averaged over the 
illustrative solid lines. 



current in the 2-direction (Fig. 4(b)). Once the sample 
is placed in a large magnetic field (/? 1), a substan- 
tial current flow develops in the z-dircction. As shown in 
Figs. 4(c) and (d), the current enters the sample along 
the top surface, flows along the z-dircction in the middle 
of the sample, and then exits along the bottom surface. 
Note that the current is parallel to the equipotcntials at 
large fields, as expected, since the angle that the cur- 
rent makes with the electric field is given by arctan(/3). 
The current pattern in Fig. 4(d) is reminiscent of that in 
a rectangular homogeneous material with perfectly con- 
ducting electrodes on the left and right boundaries [24]. 
There, the current distortions were induced by a mis- 



match of the Hall resistivity between the material and 
the electrodes. In this case, the mismatch is generated 
from the change in the angle that the current spans with 
the magnetic field, rather than arising from a change in 
the properties of the charge carriers themselves. 

Experimentally, we characterize the transport in the 
z-direction of the silver chalcogenides by measuring the 
longitudinal-transverse coupling R^z, i-e. the voltage 
drop across the sample thickness divided by the total cur- 
rent entering the sample. We plot in the main panel of 
Fig. 5(a) the R^z tensor component for an n-typc Ag2Se 
sample as a function of dimensionless field /3 = lot at T = 
110 K. Here, a six-probe configuration was used to com- 
pare the voltage drop at two different probe locations, 
and allow for later correction of possible lead misalign- 
ments [R^H) - i?..(0)^^^). 

The observed finite R^z provides clear evidence of the 
predicted current distortion in the non-stoichiometric sil- 
ver chalcogenides system. Moreover, except for obvious 
changes in the magnitude, the functional dependence on 
magnetic field of R^z appears to be rather insensitive to 
the location of voltage leads. Fig. 5(b) provides the cal- 
culated Rxz components of the percolation pattern con- 
structed from 12 X 16 X 30 elements, where the averaged 
potential on each side of the sample was determined by 
sampling along the illustrated solid lines (Fig. 5(b), in- 
set). The simulated results show good consistency with 
the experimental data, as demonstrated by an increas- 
ing Rxz with increasing H without saturation up to 14 
T (lot ^ 14). While both experiment and simulation 
have a similar dependence on probe location, note that 
the theoretical R^z must decrease as we approach the 
boundaries because the ideal electrodes do not support a 
voltage drop. 

The non-zero R^z observed in experiment and simula- 
tion indicates that the inhomogeneities are macroscopic, 
with a length scale that is appreciable compared to the 
sample dimensions. To see this, it is important to stress 
that the sample is not invariant with respect to a 180- 
dcgree rotation about the field axis, indicating that the 
current paths favor a given direction in the y-z plane. 
Such symmetry breaking is not expected to occur in an 
isotropic system that is effectively infinite, but finite sized 
systems, such as our percolation pattern, can exhibit 
a particular chirality. However, not all the salient fea- 
tures of experiment are captured by our simple model. 
The theoretical R^z is purely anti-symmetric with re- 
spect to magnetic field, while the experimental Rxz con- 
tains a substantial symmetric part. This is not surpris- 
ing given that the inhomogeneities in our model are only 
two-dimensional. To generate a symmetric dependence of 
Rxz on field, we require current fiow along the z-direction 
in the absence of magnetic fields, and therefore three- 
dimensional inhomogeneities, so that the z-component 
of the current is not completely controlled by the field 
direction. Nonetheless, our simple model serves to il- 
lustrate the principle of field-induced current distortions 
across the sample thickness. 
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C. Longitudinal Magnetoresistance 



The longitudinal MR provides the most striking man- 
ifestation of sample inhomogeneities. As the distorted 
current now lies across the thickness of the sample, bend- 
ing the cquipotcntial lines (Fig. 5(a), inset), the voltage 
drop across the length of the sample will be reduced sig- 
nificantly and a negative MR can emerge. We plot in 
Fig. 6 the normalized longitudinal MR of n-type Ag2Te 
as a function of temperature and sample thickness, and 
see a pronounced negative MR in the bulk sample (400 
fim) for T > 50 K. The maximum in the negative MR at 
T = 110 K corresponds to the crossover from intrinsic to 
extrinsic semiconducting behavior, where the distortions 
of current become the strongest as the fluctuations in mo- 
bility can veer from positive to negative. However, this 
effect is diminished for thinner samples, such that there is 
no negative MR observed up to = 14 T for thickness 20 
yiirn. This demonstrates not only that current flow across 
the sample thickness is crucial for the appearance of neg- 
ative MR, but that the critical thickness of 20 fim gives a 
measure of the inhomogeneity length scale. The effects of 
spatial randomness are suppressed completely if the size 
of the system becomes comparable to the characteristic 
inhomogeneity length scale in one dimension. 

We now compare the experimental results to simula- 
tions. For the two-dimensional percolation pattern in 
Fig. 7(a) with pi/p2 = 10,000, we know from Fig. 4(c) 
that there is a flow of current away from the surface at 
high magnetic fields. Hence a negative MR would be ob- 
tained from a four-probe measurement where the voltage 
drop is measured by placing probes on the top surface 
of the sample. However, the induced current distortions 
also increase the actual resistance of the sample, and we 
find that this effect dominates, resulting in the same pos- 
itive linear magnetoresistance as would be obtained from 
a two-probe measurement (Fig. 7(b)). This is consistent 
with the experiments, where the magnetoresistance be- 
comes both positive and quasi-linear when we reduce the 
inhomogeneities in the z-direction by thinning the sam- 
ple. Similar to the case of transverse MR, the linearity 
of the longitudinal MR is derived from the Hall resis- 
tance contributed by distorted current paths, namely by 
the currents in the z-direction of Fig. 4(d). The magni- 
tude of the longitudinal MR, however, is reduced when 
we decrease pi/ P2, as the magnetoresistance always sat- 
urates at a field scale set by pxj P2- As expected, when 
we take the limit ^ in Fig. 2(d), we find that 
Ai?/i? (X Lz, so that there is no magnetoresistance in 
a truly two-dimensional system. 

In order to obtain a negative magnetoresistance, one 
requires the current to be depleted from the surface faster 
than the increase in resistance. Since experiments im- 
ply that this scenario is possible when there are inhomo- 
geneities along the z-direction, we consider a slight modi- 
fication of our toy model to include three-dimensional dis- 
order, as depicted in Fig. 7(c). We perform a four-probe 
measurement by averaging voltages along the solid lines 
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FIG. 6: Normalized longitudinal magnetoresistance vs. sam- 
ple thickness of n-type Ag2Te. The negative component of the 
magnetoresistance diminishes in thinner crystals, and the in- 
homogeneity length scale is determined by the thickness where 
the negative MR eventually disappears (20 /im). 



in the j/-direction, and then take the component of the 
voltage difference that is symmetric with respect to the 
magnetic field. Referring to Fig. 7(d), we find a substan- 
tial negative magnetoresistance when p\l pi ~ 10,000, 
confirming that three-dimensional inhomogeneities are 
responsible for the negative longitudinal MR observed 
in the silver chalcogenides. It is of interest to note that 
there is always an upturn in the negative magnetoresis- 
tance at high fields for the simple models we considered, 
but it is unclear whether this is generic for the infinite 
system. 



V. CONCLUSIONS 

In this paper, we have compared directly the gal- 
vanomagnetic transport experiments on Agj^^Se and 
Ag2-|-5Te with theoretical simulations of inhomogeneous 
conductors in both two and three dimensions. By fit- 
ting the linear transverse MR with the PL 2D random 
resistor network model, we have demonstrated that the 
silver chalcogenides contain roughly equal proportions of 
electron- and hole-like regions, and that the field scale for 
the onset of linearity is set by the width of the mobility 
distribution rather than the average mobility. Extending 
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the model to three dimensions, we find that simple real- 
izations of macroscopic disorder are able to reproduce the 
observed finite longitudinal-transverse coupling Rxz and 
the negative longitudinal MR in the silver chalcogenides. 
Moreover, visualizations of the current distortions across 
the thickness dimension illustrate the depletion of current 
at the sample surface that is crucial for the phenomenon. 
Future work is required to extend the theoretical results 
to the infinite system and to characterize precisely the 
types of disorder which give rise to the biggest magne- 
toresistive response. 
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